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Abstract 

We consider a simple Maier-Saupe statistical model with the in- 
clusion of disorder degrees of freedom to mimic the phase diagram 
of a mixture of rod-like and disc-like molecules. A quenched distri- 
bution of shapes leads to the existence of a stable biaxial nematic 
phase, in qualitative agreement with experimental findings for some 
ternary lyotropic liquid mixtures. An annealed distribution, however, 
which is more adequate to liquid mixtures, precludes the stability of 
this biaxial phase. We then use a two-temperature formalism, and as- 
sume a separation of relaxation times, to show that a partial degree of 
annealing is already sufficient to stabilize a biaxial nematic structure. 

Quenched and annealed degrees of freedom of statistical systems are 
known to produce phase diagrams with a number of distinct features pQ. 
The ferromagnetic site-diluted Ising model provides an example of a contin- 
uous transition, in the quenched case, which turns into a first-order boundary 
beyond a certain tricritical point, if we consider thermalized site dilution [2]. 
Disordered degrees of freedom in solid compounds, as random magnets and 
spin-glasses, are examples of quenched disorder, which lead to well-known 
problems related to averages of sets of disordered free energies. In liquid 
systems, however, relaxation times are shorter, and the simpler problems of 
annealed disorder are more relevant from the physical perspective. In this 
paper, we show that distinctions between quenched and annealed degrees 



of freedom are particularly relevant in statistical models of mixtures, which 
have been used to account for the elusive biaxial nematic phase of liquid 
crystals [3]. 

Uniaxial nematic phases, with the director along a single axis, have been 
fully characterized in the phase diagrams of a large number of thermotropic 
as well as lyotropic liquid crystals [2] [I]. The existence of a biaxial nematic 
structure, however, has been subjected to some debate [5]. A biaxial ne- 
matic phase was predicted by calculations for different lattice models [6], 
and has been found by Yu and Saupe [7J in the phase diagram of a ternary 
lyotropic mixture. In the ordered region of this mixture, in a phase dia- 
gram in terms of temperature and concentration of one of the compounds, 
there is a nematic biaxial structure bounded by two distinct uniaxial nematic 
structures, loosely associated with either prolate (cylinder-like) or oblate 
(disc-like) molecular aggregates. Transitions between the ordered uniaxial 
nematic and the disordered phases are discontinuous, in agreement with the 
Maier-Saupe approach, and transitions between uniaxial and biaxial struc- 
tures are continuous, with a critical line that is supposed to end at a Landau 
multicritical point. More recent experimental work has indicated the exis- 
tence of nematic biaxial structures in a certain number of new thermotropic 
compounds, formed by anisotropic banana-shaped molecules [5]. 

Although earlier theoretical calculations indicated the existence of a sta- 
ble biaxial nematic structure [5], which is also supported by a Landau-de 
Gennes expansion [8J, mean- field calculations by Palffy-Muhoray and col- 
laborators [9], using the Maier-Saupe interactions to consider a mixture of 
cylinders and discs, precluded the stability of an intermediate biaxial struc- 
ture, except under some special circumstances. A few years ago this problem 
was reanalyzed, in terms of a schematic discrete version of the Maier-Saupe 
model, in a paper by Henriques and Henriques [10J, who pointed out the 
existence and stability of a biaxial nematic phase, bordered by two critical 
lines meeting at a Landau multicritical point [SJ, in close contact with the 
experimental findings of Yu and Saupe [7J. The calculations of Henriques 
and Henriques, however, which can be carried out for any distribution of 
molecular shapes, implicitly assumed a quenched polymorphism, which may 
not be adequate for these liquid crystalline systems. 

Given the scarcity of experimental data on these biaxial phases, the appar- 
ently conflicting theoretical results [6] |TTJ , and a few recent and not so conclu- 
sive simulations [T2][l3], there is still room for revisiting the statistical prob- 
lem of a mixture of cylinders and discs. Along the lines of the work of Hen- 



2 



riques and Henriques [TU], we then perform standard statistical-mechanics 
calculations for a simple discrete version of the Maier-Saupe model, which 
we call the basic model, with the inclusion of a binary distribution of shapes 
to mimic a mixture of prolate and oblate molecules (cylinders and discs). 
We draw clear distinctions between quenched and annealed distributions of 
shapes. In the quenched case, we show the existence of a stable biaxial phase 
in the temperature-concentration phase diagram, in qualitative agreement 
with experiments, and in accordance with the work of Henriques and Hen- 
riques. In the annealed case, however, the same model system leads to a 
unstable biaxial phase, in connection with the earlier calculations of Palffy- 
Muhoray and collaborators, which provides an example of the distinctions 
between the effects of quenched and annealed disorder degrees of freedom. 

Although the results for the quenched distribution seem to be in agree- 
ment with the experimental phase diagrams, relaxation times in both liquid 
and liquid crystalline systems are relatively short, and the physics of dis- 
ordered couplings in these liquid systems should be better represented by 
thermalized variables. Even in the annealed situation, however, we should 
make an attempt to account for a certain degree of separation between re- 
laxation times. We then resorted to a two-temperature formalism [13], with 
orientational degrees of freedom and disorder variables coupled to different 
heat reservoirs, at temperatures T and T\ = nT. In this application of the 
formalism, we show that a small difference of temperatures is already enough 
to change the phase diagram of the fully annealed case, and to stabilize the 
biaxial nematic structure, which probably provides an explanation for the 
elusive nature of this phase behavior. 

In analogy to the Curie- Weiss model of ferromagnetism, the Maier-Saupe 
model is given by the energy 

n = ~ £ £ wrsf, (i) 

l<i<j<N (i,v=l,2,3 

where A is a coupling constant, the orientational degrees of freedom are given 
by 

S? v = \ (3« - <W . (2) 

where | lti\ = 1 for i = 1, 2, ... N, and {Aj} is a set of (disordered) couplings 
representing either prolate (Aj = +1) or oblate (Aj = — 1) molecular groups. 
This model is further simplified if we suppose that the director ¥j assumes 
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six values only, along the Cartesian axes, iti = (±1,0,0), (0,±1,0), and 
(0, 0, ±1), according to an early suggestion of Zwanzig. 

In the quenched case, {Aj} is a set of independent, identical, and identi- 
cally distributed random variables, associated with a probability distribution 
p(Xi). For a given configuration {Aj}, the partition function is written as 



Z q {K} 



ex p 



N ^ 



E 



l<i<j<N n,v=l,2,3 



(3) 



where j3 is the inverse of a dimensionless temperature, and we choose A — 1. 
We now use standard Gaussian identities to write the partition function as 
an integration over the independent terms of a symmetric tensor Q. In the 
thermodynamic limit, taking into account the self- averaging properties of 
this problem, we can write 



= J [dQ] exp [- 



■PNf< 



<?J i 



where 



(4) 



(5) 



The equations of minima, dg/dQ^ = 0, lead to the usual traceless ten- 
sor order parameter of the Maier-Saupe approach. We then introduce the 
standard parametrization Q xx = — (S + r))/2, Q yy = —(S — r))/2, and 
Qzz = S, and consider a double-delta probability distribution, p (A«) = 
c8 (Aj — 1) + (1 — c) 5 (Aj + 1), where the parameter c represents the con- 
centration of prolate (A = +1) molecules. We then have 



/, = i(3S* + In 2-| In 



2 exp 
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(6) 



which leads to the phase diagram of Fig. 1, with thermodynamically stable 
uniaxial (S ^ 0; i] = 0) and biaxial (S ^ 0; i] ^ 0) nematic structures, and 
two critical lines meeting at the Landau multicritical point, (3 = 4/3 and 
c = 1/2. 
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Figure 1: Temperature-concentration phase diagram for a quenched distri- 
bution of shapes. First-order transitions between the isotropic region and 
the uniaxial nematic structures (N c and N^) are indicated by the solid line. 
Dashed lines indicate second-order transitions. 



In order to make contact with the standard Landau-de Gennes phe- 
nomenology, we expand the free energy f q in the neighborhood of the Lan- 
dau multicritical point, in terms of the quadratic and cubic, I2 = Tr Q 2 and 
^3 = Tr Q 3 , invariants of the tensor order parameter, 



According to published analyses [8] of this phenomenological expansion, a 
positive coefficient of if is indeed enough to stabilize the biaxial phase. 

We now turn to the annealed case, which is associated with the canonical 
partition function 



where the sum over over {Aj} is restricted by the fixed concentrations of the 
molecular types. As A = +1 corresponds to a prolate and A = —1 to an 
oblate molecule, we have N p — N a — J2f=i an d N Q = N — N p , where N p 
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(N ) is the number of prolate (oblate) molecules, and N is the total number 
of molecules. We then introduce a chemical potential and change to the 
grand-canonical ensemble, 
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For the same basic model, with six choices of the unit vectors, we have 



(9) 



[dQ] exp [-0N<l> a ] , 



(10) 



with the grand potential per particle 



2 

MM 



2" 



i.„2- 



i in ex p - \ X $ E ^w* j E exp fypQw) 



A=±l 

which can also be written as 



(11) 



= i(35 2 + ^-^-|ln2- 

2 cosh ( \f*S cosh (- 2 ;,, | ;.s j +cosh Q ;// • | is 



(12) 



Given the thermodynamic fields T and we show that a biaxial solution 
exists but that it is thermodynamically unstable, in agreement with the cal- 
culations of Palffy-Muhoray and collaborators. In the T — \i phase diagram 
there is just a first-order boundary, at /i — 0, for T — l/j3 < 3/4, between 
two distinct uniaxial structures. In this fully annealed situation, we can also 
draw the corresponding temperature-concentration phase diagram, as shown 
in Fig. 2, in which we sketch the tie lines of coexistence of the two uniaxial 
nematic phases. 



6 



0.8 



0.6 

T 

0.4 



0.2 







III 

ISOTROPIC 


0.7B 






T 












(1.77 










— N 


















III 



0.2 0.4 0.6 0.8 

c 



Figure 2: Temperature-concentration phase diagram in the fully annealead 
case. The horizontal tie lines indicate the coexistence between two distinct 
uniaxial nematic phases. In the inset, we show an amplification of a section 
of the (narrow) coexistence region between uniaxial nematic and isotropic 
phases. There is no stable biaxial nematic structure. 



In order to make contact with the Landau-de Gennes phenomenology, we 
can also write an expansion in the neighborhood of the Landau multicritical 
point, (3 = 4/3 and /i = 0, 



--lnl2 + - 
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(13) 



In contrast to Eq. (CO), the coefficient of Ig is negative, which now precludes 
the stability of a biaxial nematic structure [8]. 

If we use standard techniques of statistical physics, it is straightforward 
to show that the Helmholtz free energy /„ of the annealed version of the 
basic Maier-Saupe model differs from the corresponding free energy f q of the 
quenched version by a term corresponding to the entropy of mixing, 



/a = / 9 + ^[clnc+(l 



c) In (1 - c)] . 



(14) 
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This is indeed a quite general result, which we have been able to show for 
a class of mean-field self- averaging disorder variables, with the insertion of 
the appropriate form of the entropy of mixing. It is interesting to point 
out a resemblance with a derivation by Mazo of a quenched Helmholtz free 
energy that includes, in addition to the expectation of the logarithm of the 
partition function, a usually inaccessible term of entropy [13]. We remark 
that quenched situations are far from true thermodynamic equilibrium, and 
that in the quenched case we do not have access to the entropy of mixing. Of 
course, the thermodynamic analysis of the Helmholtz free energy /„, which 
should be a properly convex function of the density c, leads to the same 
results of the analysis of the grand potential a , written in terms of the 
thermodynamic fields, T and \i. 

We now search for a stable biaxial phase in a situation of partial annealing, 
which may be represented by the introduction of two heat baths, at differ- 
ent temperatures, associated with the relaxation times of the orientational 
(quicker) and disorder (slower) degrees of freedom [14J. Given a configuration 
A of the slower disorder variables, we can schematically write the probability 
of occurrence of a configuration a of the orientational variables, 

P(<7|A) = -^exp[-/W(<7,A)], (15) 
No- 
where T = 1/(3 is the temperature of a heat bath, and 

Z a = Z a (A) = £ exp [-pH (a, A)] . (16) 

a 

The time evolution of Aj is given by a Langevin equation, 

f)\ rftt 

reT—'^-fsSw- (17) 

where z (t) is a multiplier associated with the chemical potential, and 

M*)»fe(0> = 2TT x 8 ij 5{t'-t') 1 (18) 

where we have introduced the temperature T\ of a second heat bath. With 
the assumption of quick and slow time scales, it is reasonable to replace 
dTi/dXi by its average value, 

«"=,/»»\ = ^, (19) 
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where 

Heft = U eff (A) = -k B T\n Tr a exp [-0H (a, A)] . (20) 

We then assume that the probability of a configuration {Aj} is given by the 
grand-canonical expression 



E(P x ,P,N,fi) 



(21) 



where 



(Px,P,N,fi)= I [dAK^exp 
J I M 



-PU (v, A) + ^ ( £A, + N 



(22) 

with n = T/T\, which resembles the number of replicas in spin-glass prob- 
lems. The application of this formalism to the basic Maier-Saupe model 
leads to a grand potential 0, with the following expansion about the Landau 
multicritical point, (3 = 4/3 and \i — 0, 



3 3 1 
— In 6 ln2 + - 
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(23) 



For n = 1, we regain the results for the fully annealed situation, given by 
Eq. (fl3l) . From the coefficient of if, we immediately see that the biaxial 
nematic phase is stable for n < 9/10, which indicates that a relatively weak 
departure from the fully annealed situation is already enough to produce a 
stable biaxial structure. 

In conclusion, we use a basic Maier-Saupe model of a binary mixture of 
cylinders and discs to investigate the stability of a biaxial nematic phase. 
In the quenched case, the biaxial phase is thermodynamically stable. In the 
fully annealed situation, it becomes unstable, but in a situation of partial an- 
nealing, represented by couplings to heat reservoirs at different temperatures, 
we may recover thermodynamical stability. 
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